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Abstract — In this paper, resource allocation for multiple- 
input multiple-output orthogonal frequency division multiplexing 
(MIMO-OFDM) downlink networks with large numbers of base 
station antennas is studied. Assuming perfect channel state 
information at the transmitter, the resource allocation algorithm 
design is modeled as a non-convex optimization problem which 
takes into account the joint power consumption of the power 
amplifiers, antenna unit, and signal processing circuit unit. 
Subsequently, by exploiting the law of large numbers and dual de- 
composition, an efficient suboptimal iterative resource allocation 
algorithm is proposed for maximization of the system capacity 
(bit/s). In particular, closed-form power allocation and antenna 
allocation policies are derived in each iteration. Simulation results 
illustrate that the proposed iterative resource allocation algorithm 
achieves a close-to-optimal performance in a small number of 
iterations and unveil a trade-off between system capacity and 
the number of activated antennas: Activating all antennas may 
not be a good solution for system capacity maximization when a 
system with a per antenna power cost is considered. 

I. Introduction 

The demand for ubiquitous service coverage and high speed 
communications has been growing rapidly over the last decade. 
Multiple-input multiple-output orthogonal frequency division 
multiplexing (MIMO-OFDM) technology is considered as a 
viable solution for addressing this issue, as it provides extra 
degrees of freedom for resource allocation. Specifically, the 
ergodic capacity of a point-to-point MIMO fading channel 
increases practically linearly with the minimum of the number 
of transmitter and receiver antennas [1], [2]. Theoretically, 
the system performance can be unlimited if the number of 
antennas at both the transmitter and the receiver is increased. 
However, the complexity of MIMO receivers limits the per- 
formance gains that can be achieved in practical systems, 
especially for small handheld devices. In order to facilitate the 
MIMO technology in practice, an alternative form of MIMO 
has been proposed where a base station (BS) with a large 
number of antennas serves one single antenna user [3] -[8]. 
This new paradigm for multiple antenna systems is known as 
massive MIMO which does not only shift the signal processing 
burden from the receivers to the BS, but also provides a 
promising system performance. As a result, massive MIMO 
system design has recently drawn much attention from both 
industry and academia [3]-[8]. 
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In [3], high throughput gains were shown to be achievable in 
both uplink and downlink for a time division duplex (TDD) 
system with massive MIMO. In [4], the authors derived an 
achievable downlink data rate for a single-user massive MIMO 
channel under the constraint of a constant signal envelope 
per antenna. On the other hand, different receiver signal 
detection algorithms were studied and channel measurements 
were carried out for facilitating the implementation of massive 
MIMO. In [5] and [6], low complexity likelihood ascent search 
algorithms were proposed for signal detection in large-scale 
MIMO multicarrier and single carrier systems, receptively. In 
[7] and [8], channel measurements for massive MIMO systems 
were reported which revealed the actual potential system 
performance of large-scale antenna systems. A substantial 
capacity gain (bit/s/Hz) was observed with massive MIMO 
compared to single antenna systems in all studies [3]-[8]. Yet, 
the advantages of massive MIMO do not come for free and 
the number of antennas cannot be unlimited. In practice, each 
extra antenna imposes an additional power cost on the power 
budget due to the associated electronic circuitries [9]. Besides, 
both the power amplifiers (PAs) and the antenna circuits are 
sharing the same finite capacity power source at the BS. 
This joint power consumption model has not been taken into 
account in the literature so far. Therefore, the results in [3]- 
[8] which are valid for unlimited numbers of antennas and 
unlimited power supply, may no longer be applicable for the 
case of finite numbers of antennas and limited power supply. 

In this paper, we address the above issues. In Section II, 
we introduce the adopted MIMO-OFDM channel. In Section 
III, we formulate the resource allocation algorithm design as 
an optimization problem and maximize the system capacity 
for communication in MIMO-OFDM systems. To this end, an 
suboptimal iterative algorithm is proposed in Section IV. In 
Section V, we show that the proposed suboptimal algorithm 
does not only have a fast convergence, but also achieves a 
close-to-optimal performance. In Section VI, we conclude with 
a brief summary of our results. 



II. MIMO-OFDM System Model 

In this section, after introducing the notation used in this 
paper, we present the adopted system model. 




(BS) 

Fig. 1 . Illustration of a MIMO-OFDM downlink network and power supply 
model. 

A. Notation 

A complex Gaussian random variable with mean /i and 
variance a 2 is denoted by CJ\f(^i,a 2 ). = niax{0,a;}. 

[x] , = a, if x > a, [x] , = x, if b < x < a, [x] , = if 6 > 

x. [x] + = max{0,x}. |-| and ||-|| denote the absolute value of 
a complex-valued scalar and the Euclidean norm of a vector, 
respectively. <C NxAI is the space of all N X M matrices with 
complex entries and [ ] T denotes the transpose operation. 

B. MIMO-OFDM Downlink Channel Model 

We consider a MIMO-OFDM system which consists of a 
BS with multiple antennas and a mobile user equipped with a 
single antenna, cf. Figure 1 . The total bandwidth of the system 
is B Hertz and there are np subcarriers. Each subcarrier has 
a bandwidth W = B/tif Hertz. We assume a TDD system 
and the downlink channel gains can be accurately obtained by 
measuring the uplink channel based on channel reciprocity. 
The channel impulse response is assumed to be time invariant 
(slow fading). The downlink received symbol at the user on 
subcarrier i G {1, . . . , n F } is given by 

Hi = \/Ptlgh.JfiXi + Zi, (1) 

where x%, Pi, and ii <E C Nr i xl are the transmitted symbol, 
transmitted power, and precoding vector for the link from the 
BS to the user on subcarrier i, respectively. Nt ( is the number 
of active antennas allocated on subcarrier i for transmission, 
hj 6 C Nt < x 1 is the vector of multipath fading coefficients 
between the BS and the user. The elements in are assumed 
to be independent and identically distributed (i.i.d.). I and g 
represent the path loss and shadowing between the BS and the 
user, respectively. Zi is additive white Gaussian noise (AWGN) 
on subcarrier i with C7V(0, NqW), where No is the power 
spectral density. 

III. Resource Allocation Algorithm Design 

In this section, we introduce the adopted system perfor- 
mance metric and formulate the corresponding resource allo- 
cation problem. 



A. Instantaneous Channel Capacity 

In this subsection, we define the adopted system perfor- 
mance measure. Given perfect channel state information (CSI) 
at the receiver, the channel capacity between the BS and the 
user on subcarrier i with channel bandwidth W = B/rip is 
given by 

/ \ P lo\h T i I 2 

d = Wlog 2 (l + I\) and - * ™ ^ , (2) 

where I\ is the received signal-to-noise ratio (SNR) at the user 
on subcarrier i. 

The aggregate system capacity is defined as a weighted sum 
of the number of bits per second successfully delivered to the 
mobile users (bits-per-second) and is given by 

U(T,A) = Y, C » (3) 

where V and A are the power allocation policy and the antenna 
allocation policy, respectively. 

B. Power Consumption Model 

We model the power dissipation in the system as the sum 
of three terms [10] which can be expressed as 

U TP {V,A) = U (A) + U (V) + P where (4) 
U(A) = max {N^ } x Pa c and 

i 

Circuit power consumption of all antennas at the BS 

j i.p- 

U(V) = YleP .(5) 

Power consumption of all power amplifiers at the BS 

Pa c i s m e constant circuit power consumption per antenna, 
which includes the power dissipation of the transmit filter, 
mixer, frequency synthesizer, and digital-to-analog converter 
and is independent of the actual transmitted power. In the 
considered system, we assume that there is a maximum and a 
minimum number of active antennas, i.e., N^^ and iV m i n , at 
the BS. However, we do not necessarily activate the maximum 
number of antennas for the sake of efficient communication 
and the optimal number of active antennas will be found in the 
next section based on optimization. The physical meaning of 
the term maxfJVyJ in (4) is that if an antenna is activated, 

i 

it consumes power even if it is used only by some of the 
subcarriers. In other words, the first term in (4) represents 
the total power consumed by the activated antennas. The 
second term in (4) denotes the total power consumption in 
the radio frequency (RF) PAs of the BS. e > 1 is a constant 
which accounts for the inefficiency in the power amplifier and 
the power efficiency is defined as |. Pq is the basic signal 
processing power consumption which is independent of the 
number of transmit antennas. 



C. Optimization Problem Formulation 

The optimal resource allocation policies (V*,A*) can be 
obtained by solving 



max U(V,A) 



s.t. CI: 



(6) 



C2: U T p{V,A) < Ppg, C3: P t > 0, Vi, 

C4: 7V Tj G {iV min ,iV min + l,^ min + 2,...,7V max },Vi. 

F max in CI is the maximum transmit power allowance which 
puts a limit on the transmit spectrum mask to control the 
amount of out-of-cell interference in the downlink. C2 is 
imposed to guarantee that the total power consumption of the 
system is less than the maximum power supply from the power 
grid, Ppo, cf. Figure 1. C3 are the non-negative constraints 
on the power allocation variables. C4 are the combinatorial 
constraints on the number of activated antennas. 

Note that the above optimization problem formulation can 
be extended to the case of energy efficiency maximization 
with imperfect CSI and multiple users, as shown in the related 
journal paper [11]. 

IV. Solution of the Optimization Problem 

The optimization problem in (6) is a mixed combinatorial 
and non-convex optimization problem. To obtain an optimal 
solution, an exhaustive search is needed which is computa- 
tional infeasible for ,np 1. To strike a balance between 
system performance and solution tractability, we transform the 
problem transformation in the next section. 

A. Problem Transformation 

We assume that the BS chooses the beamforming vector 
to be the eigenvector corresponding to the maximum eigen- 
value of hihj, i.e., ii = The adopted beamforming 
scheme is known as maximum ratio transmission (MRT). 
Then, we introduce the following proposition by exploiting 
the properties of large numbers of antennas. 

Proposition 1: The asymptotic channel capacity between 
the BS and the user on subcarrier i for — > oo with MRT 
can be approximated by 



Ci 



where Ti 



wn og2 (i + rV) 
PjlgN Tl 
WN ' 



(t>) 



W\og 2 (iV) 



(7) 
(8) 



Corollary 1: In the limiting case of large numbers of anten- 
nas with MRT, the resource allocation algorithm is a chunk- 
based allocation policy. Specifically, 



N*. = N T , Vi, 
P.* = P, Vi, 



(9) 
(10) 



where P* and iVJ. denote the optimal power allocation and 
antenna activation solutions in subcarrier i, respectively. In 
other words, the BS has no incentive to allocate more resources 
to a particular subcarrier, since all the subcarriers have the 
same channel gain for 2Vr 4 — > oo. 

On the other hand, we need to handle the combinatorial 
constraint on the antenna activation. We relax constraint C4 
such that Nt ( can be a real value between N m i B and N max 
instead of an integer value, i.e., N m - m < < N max ., Vi. 
This yields the following relaxed problem: 



max 
v 



ax > W logo 

,X f-f B2 V WN J 



(ID 



s.t. CI: J2 p i < p max, C2: U T p{V, A) < P PG , 

i=l 

C3: > 0, Vi, k, C4: iV rain < N Ti < N max , Vi. 

Now, we focus on the resource allocation algorithm design 
for the relaxed problem. Note that the system performance of 
the relaxed problem will serve as an upper bound of (6). It 
can be shown that the relaxed optimization problem is jointly 
concave with respect to (w.r.t.) Pi and N^, cf. Appendix. 
Besides, it can be shown that the relaxed problem satisfies 
Slater's constraint qualification. As a result, for the relaxed 
problem, solving the dual problem is equivalent to solving the 
primal which is known as strong duality [12]. 

B. Dual Problem Formulation 

In this subsection, we solve the power allocation and 
antenna allocation optimization problem by solving its dual. 
For this purpose, we first need the Lagrangian function of the 
primal problem. Upon rearranging terms, the Lagrangian can 
be written as 



C(X,(3,V,A) 



i=i 



WN 



a) and (b) are due to the law of large numbers and —/3^Nt X Pa + ePi + Pq — Ppc 



oo, i.e., limjv T .^oo -j^- 1 = 1 and log 2 (l + x) w log 2 (x) for 
x ^> 1, respectively 1 . 

It can be observed from (7) that all the subcarriers have 
the same asymptotic channel gain due to channel hardening 
[1]. Then, we have the following corollary for the considered 
resource allocation algorithm design. 



(12) 



'in practice, the value of N n 
numbers holds. 



is selected such that the law of large 



We note that max{Nr i } has been replaced by in (12) due 

i 

to Corollary 1 . A > is the Lagrange multiplier corresponding 
to the maximum transmit power allowance CI. /3 is the 
Lagrange multiplier accounting for the power usage from the 
power grid. On the other hand, boundary constraints C3 and 
C4 will be absorbed into the Karush-Kuhn-Tucker (KKT) con- 
ditions when deriving the optimal resource allocation policies 



Subgradient update of A.,f} according to (16), (17) 





2 nd layer: Master Problem 
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Information passing from Subproblems to Master Problem 

^ Information passing from Master Problem to Subproblms 

Fig. 2. Illustration of the dual decomposition of a large-scale problem into 
a two-layer problem. 



in the following. Thus, the dual problem is given by 



mm max 

\,/3>o v,A 



(13) 



C. Dual Decomposition 

By Lagrange dual decomposition, the dual problem is 
decomposed into two layers: the first layer consists of nj? 
subproblems with identical structure; the second layer is the 
master problem, cf. Figure 2. Then, the dual problem can be 
solved iteratively, where in each iteration the transmitter solves 
the subproblem by using the KKT conditions for a fixed set of 
Lagrange multipliers, and the master problem is solved using 
the gradient method. 

1 ) Subproblem Solution: Using standard optimization tech- 
niques and the KKT conditions, the closed-form resource 
allocation policies in subcarrier i are obtained as: 

i + 

P* = P = 



N't 



N 7 



ln(2)(A + /3e) 



W 



\n{2)P Ao p 



N„ 



Vi, 



,Vi. 



(14) 



(15) 



J N„ 



The resource allocation solution in (14) and (15) can be 
interpreted as a water-filling scheme. It can be observed 
that the BS does not always activate all the antennas for 
maximization of the system capacity which is limited by the 
total power supply Ppc via /3, cf. (15). Note that since the 
solution of the np subproblems are identical, the complexity 
in solving the np subproblems can be reduced by a factor of 
np by just focusing on one subproblem. 

2) Master Problem Solution: The dual function is differ- 
entiable and, hence, the gradient method can be used to solve 
the second layer master problem in (13) which leads to 



A(m 

f3{m 



1) 



D = 



A(m) - £i(m) x (P mi 

0(m)-&(m) 

X {PpG -N T x P Ac 



n F P) 



n F eP - P 



(16) 



(17) 



where index m > is the iteration index and £ u (m), 
u e {1,2}, are positive step sizes. Since the relaxed prob- 
lem is concave in nature, it is guaranteed that the iteration 
between master problem and subproblems converges to the 
optimal solution of (13), if the chosen step sizes satisfy the 
infinite travel condition [12], [13]. Then, the updated Lagrange 
multipliers in (16), (17) are used for solving the subproblems 
in (13) via updating the resource allocation policies. 

Note that the resource allocation solutions obtained in (14) 
and (15) are optimal w.r.t. the relaxed problem in (11) for 
N Ti > 1. 

V. Simulation Results 

In this section, we study the performance of two algorithms 
via simulation. The first algorithm achieves the upper bound 
performance of (6) by solving the relaxed problem in (11) as 
outlined in Section III-C. For the second algorithm, we apply a 
floor function |_-J to the antenna allocation solution in (15), i.e., 
= L-^r J ' which achieves a suboptimal performance w.r.t. 
the original problem formulation in (6). For the system setting, 
there are np = 128 subcarriers with carrier center frequency 
2.5 GHz and a total system bandwidth of B = 5 MHz. We as- 
sume a noise power of NqW = — 118 dBm in each subcarrier. 
The desired user is located at a distance of do = 500 m from 
the BS. Log-normal shadowing with a standard deviation of 8 
dB is assumed. The multipath fading coefficients of the BS-to- 
user link are modeled as i.i.d. Rayleigh random variables with 
unit variances. The average system capacity is obtained by 
assuming capacity-achieving codes and counting the amount 
of data which is successfully decoded by the user averaged 
over both shadowing and multipath fading. We assume a 
static circuit power consumption of Pq = 40 dBm [14]. The 
additional power dissipation incurred by each extra antenna 
for transmission is Pa c = 30 dBm [9]. On the other hand, 
we assume a power efficiency of 40% in the RF PA, i.e., 
e = jpj = 2.5. The maximum and minimum numbers of active 
antennas are set to N max = 500 and iV m i n = 10, respectively. 



A. Convergence of Iterative Algorithm 

Figure 3 illustrates the evolution of the proposed iterative 
algorithm for different values of maximum transmit power al- 
lowance, F max , and maximum power supply, Ppg- The results 
in Figure 3 were averaged over 100000 independent adaptation 
processes where each adaptation process involves different 
realizations of the shadowing and the multipath fading. It can 
be observed that the iterative algorithm converges to values 
close to the upper bound performance within 4 iterations for 
all considered cases. In other words, a close-to-optimal system 
capacity can be achieved within a few iterations on average. 

In the following case studies, we focus on the suboptimal 
algorithm and set the number of iterations in the algorithm to 
5. 

B. Average Capacity versus Maximum Power Supply 

Figure 4 illustrates the average system capacity (bit/s/Hz) 
versus the maximum power supply from the power grid, Ppa, 




Fig. 3. Average system capacity (bit/s/Hz) versus number of iterations 
with different maximum powers supplied by the power grid, Ppc, and 
different values of maximum transmit power allowance, P max - The dashed 
lines represent an upper bound on maximum average system capacity for 
different cases. 

for the proposed suboptimal algorithm. Different values of 
maximum transmit power allowance, P max , are considered. It 
can be observed that the average system capacity increases for 
increasing amount of power supplied by the power grid, since 
more power is available at the BS for resource allocation. On 
the other hand, the system performance also benefits from in- 
creasing values of P ma x- Indeed, from an optimization theory 
point of view, a larger value of P max spans a larger feasible 
solution set which results in a possibly higher objective value. 
However, there is a diminishing return in performance as P max 
increases, especially in the low power supply regime, i.e., 
Ppg < 49 dBm. This is because the system performance is 
limited by the small amount of power Ppc supplied by the 
power grid, instead of P max . 

C. Average Number of Activated Antennas and Power Con- 
sumption versus Maximum Power Supply 

Figure 5 depicts the average number of activated antennas 
versus maximum power supply, Ppg, for the proposed subop- 
timal algorithm. The number of activated antennas increases 
rapidly w.r.t. the an increasing Ppg due to a larger amount of 
available power at the BS. However, the curves in Figure 5 
have different slopes for different values of P max - In particular, 
the number of activated antennas increases with Ppc much 
faster for a small value of P max than for a large value of 
Pmax- This is because for a small value of P max , the system 
performance is always limited by the amount of power radiated 
in the RF. As a result, the BS tends to activate more antennas 
for increasing the antenna array gain. On the other hand, when 
the value of P max increases, the BS will allocate more power 
to the RF and decrease the number of activated antennas. 
Counterintuitively, the BS does not activate all antennas for 
maximizing the average system capacity, when both the PAs 
and the antennas are consuming powers from the same power 



Fig. 4. Average system capacity (bit/s/Hz) versus the maximum power 
supplied by the power grid, Ppc , for different values of maximum transmit 
power allowance, Pmax- 
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Fig. 5. Average number of activated antennas versus the maximum power 
supplied by the power grid, Ppa , for different values of maximum transmit 
power allowance, P max - 



source. 

Figure 6 depicts the power consumption ratio between the 
PAs and antennas circuits, i.e., , versus maximum power 
supply, Ppc for the proposed suboptimal algorithm. We can 
see that for a small value of P max , the BS spends more 
power on activating antennas when Ppc increases. This is 
because more antennas will be activated to enhance the system 
capacity when the radiated power in the RF reaches P max - 
This observation coincides with the results in Figure 5. On 
the other hand, when P max increases, the power consumption 
of both the PAs and the antenna activation approach the same 
value for maximizing the system capacity. Indeed, the BS does 
not have a higher preference for spending more power in the 
PAs than for activating more antennas, or vice versa, when 
the constraints on both P max and N max in (6) become less 




Fig. 6. Average power consumption ratio between PAs and antennas circuits, 
i.e., versus the maximum power supplied by the power grid, Ppa, for 

different values of maximum transmit power allowance, Pmax- 



stringent. On the contrary, the BS tries to maintain a balance 
between the power consumption in the PAs and the antenna 
circuits. As a result, activating all available antennas may not 
be a good solution for system capacity maximization, since 
the power consumptions in the antennas and PAs should be 
balanced. 

VI. Conclusion 

In this paper, we formulated the resource allocation algo- 
rithm design for MIMO-OFDM networks with large numbers 
of BS antennas as a non-convex optimization problem, where 
a joint power consumption model for the antennas, power 
amplifiers, and signal processing circuit unit was taken into 
consideration. An efficient iterative resource allocation algo- 
rithm with optimized power allocation and antenna allocation 
policies was proposed. Our simulation results did not only 
show that the proposed algorithm converges within a small 
number of iterations, but also unveiled the trade-off between 
maximum system capacity and the number of activated anten- 
nas. 



Since A r < 0, r = {1,2}, so H(G\) is a negative semi- 
definite matrix and Ci is jointly concave w.r.t. Pi and N^- 
The summation of Ci over i preserves the concavity of the 
objective function in [12]. On the other hand, constraints Cl- 
C5 in (1 1) with relaxed span a convex feasible set and thus 
the transformed problem is a concave optimization problem. 
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Appendix - Proof of Concavity of Relaxed 
Problem 

Let H(Ci) and Ai,A2 be the Hessian matrix of function 
Ci and the eigenvalues of H(Ci), respectively. The Hessian 
matrix of function Ci can be expressed as 



H(Q 



—w 

P'f ln(2) 

1 





-w 

N%. ln(2) 



(18) 



Therefore, the corresponding eigenvalues of the Hessian ma- 
trix are given by 



Ai 



-W 
7fln(2) 



and Aq 



-W 
/Vf . ln(2) ' 



(19) 



